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We study the cosmological evolution of massless single-field scalar-tensor theories of gravitation 
from the time before the onset of e^e“ annihilation and nucleosynthesis up to the present. The 
cosmological evolution together with the observational bounds on the abundances of the lightest 
elements (those mostly produced in the early universe) place constraints on the coefficients of the 
Taylor series expansion of a{ip), which specifies the coupling of the scalar field to matter and is the 
only free function in the theory. In the case when a(ip) has a minimum (i.e., when the theory evolves 
towards general relativity) these constraints translate into a stronger limit on the Post-Newtonian 
parameters 7 and /3 than any other observational test. Moreover, our bounds imply that, even at 
the epoch of annihilation and nucleosynthesis, the evolution of the universe must be very close to 
that predicted by general relativity if we do not want to over- or underproduce '^He. Thus the 
amount of scalar field contribution to gravity is very small even at such an early epoch. 


I. INTRODUCTION 

Scalar-tensor (ST) theories of gravitation have a long history. The first ones to appear were those of Jordan [||, 
Fierz [||, and Brans-Dicke |^. These are the simplest theories in the sense that they consist of just one massless scalar 
field and its coupling strength to matter is constant. Later Bergmann Q, Nordtvedt ||^, and Wagoner generalized 
the theory. In their versions the scalar field has a dynamic (scalar field dependant) coupling to matter and/or an 
arbitrary self-interaction. More recently, in a very comprehensive work, Damour and Esposito-Farese considered 
the theory with multiple scalar fields. 

One could say that scalar-tensor theories of gravity are the simplest alternative to Einstein’s general theroy of 
relativity (GR). Most recent attempts at unified models of fundamental interactions, i.e.,string theories |^, predict 
the existence of scalar partners to the tensor gravity of GR, and would have ST gravity as their low energy limit. In 
addition, these theories are important in inflationary cosmology where they provide a way of exiting the inflationary 
epoch in a non-fine-tuned way . Scalar-tensor theories provide an important framework for comparison with results 
of gravitational experiments that are going to take place in the near future, e.g Gravity Probe-B LIGOS [p] . 
Geo 0, and VIRGO 

In the present paper we study the cosmological evolution of the metric and scalar gravitational fields in ST gravity. 
We consider the theories where the scalar field is massless, i.e., no self interaction. Although it seems natural that 
without a gauge invariance the scalar would acquire a mass, certain results in string theory suggest that massless scalars 
arise naturally HU- Furthermore, it has been shown that when the scalar field is massless, the cosmological 
evolution drives the scalar field to a value that minimizes the coupling function, towards making the contribution of 
the scalar field very small. This effect makes the theory practically indistinguishable from GR at the present epoch. 

Solar system tests can put limits on the Post-Newtonian parameters 7 and /3 ||^, which in ST gravity translate 
into limits on oi and 02 : the first two coefficients of the Taylor series expansion about today of a{{p), the coupling 
function of the scalar field (p to matter. We have the relations |5|j^: 
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From analysis of the Shapiro time delay in 
have 


the Viking Mars data [|^, and from VLBI radio wave deflection |19 


we 


af<lx 10-^ (1.2) 

while Lunar Laser Ranging |2^ gives the limit 

|1 +a2|ai < 5 X 10"^. (1.3) 

These limits (the most restrictive known to the authors) imply that at the present epoch in the weak field regime 
gravity’s behavior does not differ much from that predicted by GR. We will show that the cosmological constraint 
places stronger limits on oi if 02 is positive. In addition it implies that if there is a scalar field contribution to gravity, 
it must be very small from the epochs of e“'"e“ annihilation and nucleosynthesis up to the present. The cosmological 
constraint has been analyzed less generally by Casas et. al. Kalligas et. al. [Q, and Serna and Alimi |^. A 
preliminary version of our results appears in p4[| . 

Throughout this paper we follow the convetions of MTW |^(Go = c = 1). 

II. SCALAR-TENSOR THEORIES OF GRAVITY 

We consider a theory of gravitation in which the gravitational interaction is carried by the ususal spin-2 field, the 
metric plus a spin-0 field, a scalar field ip. We work in the Einstein conformal frame of |Q,|l6|, which we call the 
spin-frame, because the spin-2 and spin-0 excitations are then decoupled from each other. The main advantage of the 
spin-frame is that, in general, calculations of the behavior of the fields (as opposed to the motion of matter) are more 
readily manageable. Another advantage of the spin-frame is that the Cauchy problem is well posed 0. 

The field equations for our theory follow from extremizing the action 

S'= J d‘^x^/^[R-2g^''d^ipdi,ip] + Sm- (2.1) 

Here g = detg^jy, R = g^'^R^y is the Ricci scalar, and Sm is the matter action, which includes all the matter fields 
and their interactions including those with the gravitational fields. The matter fields couple to the gravitational fields 
through the physical metric 

g^iu = ( 2 . 2 ) 

From the previous equation we see that the scalar field couples conformally to matter via the function A{ip). A 
priori A(ip) is a completely unspecified function of the scalar field. Solar system experiments constrain this function 
somewhat, and we will show that cosmology can give additional strong restrictions on this function. 

The three functions 


a{p) = In A((/j) 

(2.3) 

da{p) 

(2.4) 

daip) 

(2.5) 


will be very useful because they are related to observable quantities. It will be convenient to express these functions 
as their Taylor series expansion about the asymptotic cosmological value of the scalar field in the present epoch, ipo'. 


a{ip) =ai{(p- ipo) + -^(v? - H- 

a{(p) = ai+ 02(^3 - Po) - 

k {( p ) = 02-1 -, 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 


where we scale our coordinates so that a{(po) = 0. 

When we vary the action with respect to and p, we obtain the field equations 
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where 


8T^Tfii, + 2 {p,fi‘p,u - 

(2.9) 

Uip = —A^TTa{(p)T , 

(2.10) 

2 5Sm 



(2.11) 

T = 

(2.12) 

Op = -fi 

(2.13) 


are the spin-frame stress-energy tensor, its trace, and the divergence of the gradient of ^p. The Bianchi identity leads 
to the equations of motion 


,, = a{p)Tip,^. (2.14) 

In all these formulae we follow the usual convention that comas represent partial derivatives and semi-colons represent 
covariant derivatives. 

From the field equations we see that if we set p = K, with K a constant such that a{K) = 0 and 

(GFC) 

gill/ a solution to Einstein’s equations of GR, we have a solution to our field equations. It follows that all general 
relativistic solutions are also solutions in ST gravity with a appropriate choice of constant scalar field. Moreover 
since a{ip) = K locates an extremum of our coupling function a((/?). When K produces a minimum of aipp), the 
GR solutions are stable solutions of ST gravity, while if K produces a maximum or a inflection point the solutions 
are unstable. It has been demonstrated |I6| that if a{if) has a minimum or approaches one symptotycally, then the 
cosmological evolution of during the matter era drives the scalar field towards the minimum value of a, and therefore 
the universe evolves towards general relativity. 

We can express ||^,|^ the measured gravitational constant far from all matter, G, and the PPN parameters 7 and 
[3 in terms of the Taylor expansion coefficients of a{ip) 


G = l + al 

(2.15) 

-2al 

7 — 1 = - 0 

1 + 

(2.16) 

1 ala2 

^ 2(l + a2)2- 

(2.17) 


From these last formulae we see how solar system tests give us information about the coupling function A{ip)^ i.e., by 
constraining the possible values of ai and 02. Our cosmological constraint will give us information on these same two 
Taylor series coefficients which, in turn, can be translated into a restriction on 7. 


III. SCALAR-TENSOR COSMOLOGIES 

We now turn to the study of the evolution of the universe as a whole in ST gravity. We assume the universe is 
homogeneous and isotropic, so that the spacetime has a Friedmann-Robertson-Walker metric 

ds^ = -dt^ -I- R^{t) [ 2 > (3-1) 

with dU, being the element of solid angle and k = —I, 0, or 1 according to whether the universe is open, flat, or 
closed. In order to compare with observations we need to express all the physically meaurable quantities (proper time, 
expansion rate, den sity, pressure, etc.) in terms of spin-frame quantities (the ones we use for our analysis). Because 
of the relation (|2.2|) between the physical metric and the spin-frame metric, the proper time measured by ideal clocks 
is not dr = \/—ds'^; it is 


df = A{ip)\/—ds^ . (3.2) 

In particular, for comoving observers df = A{ip)dt. We shall call r the spin-frame “proper” time, but rememeber that 
it is not the time measured by ideal comoving clocks. The physical scale factor of the universe is 
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Rit) = A{ip)R{t) (3.3) 

instead of R{t), the spin-frame scale factor. The condition of homogeneity implies that our scalar field is only a 
function of the the spin-frame time coordinate, tp = ^>(1). 

The matter in the universe behaves to a large degree like a perfect fluid 

= {p + p)u^u''+ . (3.4) 


In the previous equation 


u 


M = 


dxf^ 

dr 


(3.5) 


is the spin-frame 4-velovity normalized to = — 1 for massive particles. Since dr = A{ip)dT, we have = 

A{ip)u>^ with 


dx>^ 

df 


(3.6) 


being the physical 4-velocity, i.e., the actual 4-velocity particles have. The physical stress-energy tensor is defined by 
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(3.7) 


analogous to equation (2.11) for the spin-frame stress-energy tensor. Using the relation (2.2) between the physical 
metric and the spin-frame metric, we get 




(3.8) 


which leads us to the following expressions for the spin-frame density and pressure of the matter in terms of the 
physical quantities, 


P = A'^{tp)p, p = A^{(p)p. 


(3.9) 


We remind readers that quantities with a tilde are the physical ones, which are the usual measured values. 

Substituting the metric for our isotropic and homo gene ous cosmological model and the perfect fluid spin-frame 
stress-energy tensor into our field equations (^) and (|2.10|), we obtain 


- 3^ = 47r(p -k ip) + 2((/))2 , 

®(i) 

R 

ip + i—p = -47ra((/j)(p - ip ), 


(3.10) 

(3.11) 

(3.12) 


where the overdots denote derivatives with respect to spin-frame coordinate time, t. From the Bianchi identity (2.14) 
we obtain the “conservation” equation 


d{pR^) +pd{R^) = {p — ip)R^da{(p ), 


(3.13) 


which, of course, is a consequence of the field equations. 

Just as is the case in GR, in ST gravity we need an equation of state in order to solve our field equations. In general, 
equations of state relating density and pressure are valid only for physical quantities, p = p{p), and not spin-frame 
quantities. All the cosmological epochs we consider, except the e“'"e“ era, have simple equations of stat e of the type 
p (X p. We will have a similar type of relation between the spin-frame quantities because of equation (3.9) and the 
simple form of the equation of state. The e+e“ era is handled by applying first order perturbation theory in the effect 
of the scalar field (arguing why it is not necessary to go to higher orders). To first order the spin-frame density and 
pressure are equal to the physical ones, enabling us to use the same equation of state for the spin-frame quantities. 
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IV. COSMOLOGICAL EVOLUTION AND NUCLEOSYNTHESIS 


We now turn our attention to the question of what cosmology says about ST gravity. As expected, because of 
the presence of the scalar field, the expansion rate of the universe will be different than the one predicted from GR. 
Nucleosynthesis is sensitive to the speed-up factor (the ratio of the expansion rate at nucleosynthesis to that 
predicted by GR). Thus we see that the production of primordial light elements, especially '*He, will provide us with 
limits on ST gravity. In order to obtain these constraints we will integrate semi-analytically (in perturbation theory) 
our cosmological equations from the present era back to the epoch of e“'"e“ annihilation. This last era is handled 
numerically, performing the appropriate matching at the transitions between eras. 

The next step is to establish the relation between the speed-up factor, the scalar field, and the other quantities 
of our theory. The following relation between the speed-up factor, and the conformal factor, A{ip), has been 
derived [|^ from quite general considerations. 


A{ipn) 

yjl + a‘^{(po) ’ 


(4.1) 


with ipn being the value of the scalar field at the time of nucleosynthesis and ipo its value today. Since cannot differ 
from unity by more than ~ 15% without over- or under producing ^He pq , p^ we can Taylor expand our coupling 
function about ipo and consider terms beyond the second order as unimportant for all times between the present and 
nucleosynthesis: 


a{ip) = ai{ip - ipo) + y((/5 - Vof ■ 


Combining the two previous equations with equation (2.3) we arrive at 


/ ,2 1 2 


In^n = ai(v3„ - ipo) + y (‘^n - ^o) - -^ai , 


(4.2) 


(4.3) 


to the desired order. 

The equations of motion can be integrated to yield ipn — ^o- When this quantity is substituted in our last equation 
we will have the speed-up factor as a function of the the Taylor coefficients oi and 02. Then we can use the abundances 
of primordial light elements ( functions of ^„) to obtain the allowed values of oi as a function of 02. 

We will restrict ourselves to the study of those ST cosmologies that evolve towards GR. This is equivalent to 
considering only coupling functions a{(f) that at least have a minimum or approach one asymptotically. In terms of 
the Taylor expansion coefficients of the coupling function, this implies 02 > 0. Excluded is the Brans-Dicke case [a((/3) 
linear in (/?], which is not that interesting because the deviations from GR are fixed, i.e., a does not depend on the 
scalar field. Excluded also are the cases where the coupling function has no minimum, has points of inflection and at 
most one maximum. In this case solar system constraints place us very near a point of inflection or the maximum. In 
either case the coupling function will evolve away from such a point and, in the future, the deviations from general 
relativity will be more pronounced. We will consider both fiat and open universes. For each universe we will have the 
following eras: e“'"e“ annihilation, radiation dominated, and matter dominated. For the open universe we also have a 
curvature dominated epoch. We approximate the transitions between eras as sudden. 


A. Scalar Field Perturbation Theory 


Scalar-tensor theories are very general because the coupling function a(}p) could be any function of the scalar field. 
This generality prevents one from performing any calculation unless one chooses a particular coupling function. Thus 
it seems that any calculations will be dependent on the choice among the uncountably many coupling functions. 
This appears quite unnatural especially since the coupling function a{}p) should be a fundamental property of nature 
(a constant corresponding to GR). Therefore a((p) should be derived from fundamental principles or somehow be 
obtained from experimental facts. We try to do the latter in the present work using the light element abundances. 
The function a((p) = ^ acts as a measure of the deviations from GR. When |q;| ^ 1 the scalar contribution to 


gravity is small and the theory is very close to GR, while if a > 1 there wi ll b e large dev iations from GR. The solar 
systems limits on the PPN parameters 7 and j3 mentioned in equations (ED and ( |1.3D imply that the value of a 
today, Cl, is small. So at the present time the universe has small deviations from a general relativistic behavior. This 
suggests that we can treat the effects of the scalar contribution as a first order perturbation on GR. Since we are 
very near a minimum of 0(95), ai acts as a “measure” of how far we are from that minimum, i.e., the strength of the 
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scalar contribution to gravity. For ai = 0, we have no scalar contribution (i.e., GR) while for |ai| <C 1 we have small 
scalar contributions. At present, since the effect of the scalar field is small, we expect |oi| <C 1 and the cosmological 
field today, (po, to be different by order oi from a value that minimizes a((p), i.e., a((p) = 0. Moreover, as long as we 
can apply first order perturbation theory we expect the deviations of the scalar field from its value today to also be 
proportional to ai: 


ip = ifio + aipi + O(a^). 


(4.4) 


Therefore we can consider oi as the small parameter for a first order perturbation theory expansion about ge nera l 
relati vity. Besides (p itself, a is the only other scalar field dependent quantity to appear in the field equations (2.E) 
and (2.10). Hence to do first order perturbation theory in the effect of the scalar field we need a to linear order in p: 


a(p) = Oi + a 2 (p - Po) + 0(2) = ai(l + a 2 Pi) + 0{al). (4.5) 

Since this is the derivative with respect to p of equation ( [4.2[ ), we see that we need a{p), or equivalently A(p), to 
second order. But if higher order deviations are important, we will over- or underproduce ^He. Therefore our first 
orde r perturbation theory approach is valid at least up to the time of nucleosynthesis. If we look at the field equations 
(2.9) for the metric, we see that they depend on derivatives of p squared. Since the physical metric as given by 
equation (|2.2|) differs from the spin-frame metric by quantities of second order (af), it must differ from the GR metric 
by quantities of at least second order. Moreover, the matter lagrangian, and hence, the stress-energy tensor, will be 
equal to the GR quantities to at least zeroth and first order. Therefore the metric will be unaffected to first order 
while to zeroth order we just have a GR solution, 


+ 0(al). (4.6) 

In general, the physical metric is of more interest. The first order correction in the scalar field affects the physical 
metric to second order. The second order correction to the metric affects ideal clocks and measuring sticks to the same 
order as the first order scalar field correction. Therefore we would need to find this metric correction to understand 
the full physic s of a proMem to the desired order. For the nucleosynthesis constraint we are fortunate because we 
have formula (0) or (|^) which relates the speed-up factor to the coupling function independently of the metric. 
We see that in this particular case we can neglect the back reaction of the metric in the presence of the scalar field. 
There will also be corrections to the spin-frame stress-energy tensor, but as we previously argued they will be second 
order: 


+ 0{al) , 


(4.7) 


(GFC\ 

with the general relativistic stress-energy tensor. To this order the spin-frame stress-energy tensor is equal to 

the physical one. _ 

Now we consider the ST cosmological and “conservation” equations ( I.10| -3.13) in perturbation theory. We note that 
t hese equations are not independent (they are related through the Bianchi identity). Hence we do not use equation 
( |3.10 ). For the spin-frame expansion factor, scalar field, density, and pressure we have the following relations to the 
desired order: 


R(t) = + 0(al) 

p = Po + aipi + 0(al) 
p(t) = p(^^\t) + Oial) 

Pit) =p^^^Ht) +Dial). 

When we substitute these in the cosmological equations we obtain to zeroth (GR) order: 

2 

+ 3 ;^ = 

d(pR^) + pd(R^) = 0, 

where we have dropped the (GR) superscripts. To first order scalar field equation becomes 

R 

Pi -I- + 47r(l -I- a2Pi)ip - 3p) = 0 . 



(4.8) 

(4.9) 

(4.10) 

(4.11) 


(4.12) 

(4.13) 

(4.14) 
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For later use we define 


ij} =1 + 


(4.15) 


so that xj) satisfies 


R 

xp + 3—xp + 47ra2(p — 3p)xj; = 0 . 
K 


(4.16) 


For the nucleosynthesis constraint we will need the speed-up factor. Substituting the perturbation expansion ( |4.9| ) of 
ip into the relation (O) for the speed-up factor we get 


1 ' ' M .2/ 1 , - 1 


ln^„ = ai[-^ + ipi{tn) + yV5?(i«)] = a?(-| + 


209 


(4.17) 


where is the time at nucleosynthesis, and we used equation ( 4.15 ) to get the last equality. 

Next we find solutions to the field equations in the different epochs that occur in the two types of universe we 
consider: flat and open. 


B. Curvature Dominated Epoch 


We consider an era in which the evolution of the universe is mainly driven by the curvature term, i.e., the curvture 
term is large compared to the density term. This describes to a fairly good aproximation the present epoch for a 

sufficiently open universe (fc = —1). _ 

When we neglect the density term in our zeroth order field equation ([4.12) we have 




= 0 , 


(4.18) 


which trivially integrates to 


i?(t) — t Cl , 


(4.19) 


with Cl a constant of integration. 

We must consi der ho w the matter density depends on the scale factor since we should include the matter density 
term in equation (4.16) for the scalar field because altho ugh i t is small at each time it builds up during the evolution. 
Taking p = 0 in the zeroth order conservation equation (4.13) gives us the well known simple result 


Pm — 


’ 


(4.20) 


where Mq is a constant of normalization and the subscript m on p is to specify that it is the matter density. 

To first order we use the auxiliary scalar field xp and its corresponding equation ( 4.16[ ). Using our last result the xp 
equation becomes 


xP- 


-xp- 


47102 Mo 


t + cr ' (t + ci)3 

This is a transformed Bessel’s equation [M, which has solution 


■0 = 0 . 


m = 


' 167rMoa2 

t + Ci 


Cs 


t + Ci 


-w 


116 nMoa 2 ^ 

,V t + Ci ) 


(4.21) 


(4.22) 


where the Ck are constants of integration, and J 2 {x) and Y 2 {x) are Bessel functions of the first and second kind 
respectively. _ 

It is interesting that xp in equation ( 4.22| ) does not converge to zero for large times, rather it approaches a constant. 
To first order xp = 0 corresponds to the value of ip for which a{p) = 0. Therefore the universe does not evolve towards 
a state in which the deviations from GR disappe^ but as it evolves the deviation from GR converges to a finite, 
though small, value. The attractor mechanism of [Q is not too effective because for large times the Hubble damping 
makes the scalar field “stop” before reaching the value which produces a = 0. The full attractor mechanism is present 
only for flat scalar-tensor cosmologies because those will have GR as a limit for large times. 
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C. Matter Dominated Epoch 


In a matter dominated epoch the evolution of the universe is driven by the density of the matter (described as 
a fluid of nonrelativistic particles). Hence the matter density is much larger than either the curvature term or the 
density and pressure of radiation. Such an epoch corresponds to the present for a flat universe (fc = 0), or to the 
epoch previous to the curvature do mina ted one in an open universe (fc = —1). 

The zeroth order field equation (4.12|) becomes 


3 



= Stt 


A4 


This equation is easily integrated to yield 




(4.23) 


(4.24) 


with Di a constant of integration. 

We need the scalar field evolu tion t o firs t order. Just as we did in the previo us se ction ( IVB ), we use the auxiliary 
field 'i/i. When the results S) and (|4.24|) are substituted in the ip equation (|4.16D we have 


iP- 


-iP 


202 


t + Di 3(t + Di) 


;1p = 0. 


(4.25) 


By making the variable change x = ln(t + Di) we end up with an equation for a damped harmonic oscillator: 

Ip” + Ip'+ ^a2ip = 0 , (4.26) 

where primes denote derivatives with respect to x. This equation has three types of solution (under, over, or critically- 
damped), according to the value of 02 . 

Case 1 (underdamped: 02 > |): 


ip{t) 


D 2 

\/t -f Di 


sin 



Case 2 (criticallydamped: 02 = |): 


ip{t) = 


D. 


y/t -f Di 


In 


t Di 
D 3 


Case 3 (overdamped: 02 < |): 


ip{t) 


T>2 

\/t D\ 


sinh 



(4.27) 


(4.28) 


(4.29) 


In all three cases Dk (k = 2 and 3) are constants of integration. 


D. Radiation Dominated Epoch 


Previous to the matter era, there is an epoch in which the density and pressure of relativistic particles (radiation) is 
large compared to the density of the matter (nonrelativistic fluid) as well as the curvature term. This epoch describes 

the universe fairly well prior to matter-radiation equality. _ 

Remembering that the equation of state for radiation is Pr = Pr/3^ our zeroth order conservation equation ( 4.13| ) 
becomes 


0 = d{prR^) + ^d{R^) = 


diprR^ 


(4.30) 

































which trivially gives 


Pr — 




Pr = 


No 

3i?4 


(4.31) 


with No a integration constant. The subscript r specifies that these a re th e density and pressure of radiation. 
Taking into account this last result, the expansion factor equation ([1.12|) to zeroth order is 


= Stt 




We solve this equation to find 


, , / 32ttNo \ ^ , ,1 

Rit)=i ^—^—j {t + Bi)2 


(4.32) 


(4.33) 


To first order we must solve equation (4.16) for tjj. The last term in this equation is proportional to p — 3p, which 
vanishes if we consider radiation. On the other hand the matter density {pm = also contributes to the total 
density. This is a small correction at any particular time in the radiation epoch, but it builds up as we integrate the 
■0 equation over the whole radiation era. Therefore we must include the effect of the matter density in the scalar field 
equation. We find 


0 ■ 




a 2 X 


2 {t + Bi) 


-0 = 0 , 


where we have introduced 


X = 


3Mo f 2 t:N o\^ 
4W 1“ 


This is a transformed Bessel’s equation 


which has solution 


0(0 = 


B. 


{t + Bi)2 


-Ji[4(a2x)2(t+ Bi)' 


B. 


{t + Bi )2 


rTi[4(a2x)K^ + ^i)' 


(4.34) 


(4.35) 


(4.36) 


where Ji{x) is a Bessel function of the first kind and yi{x) is a Bessel function of the second kind. In addition 
Bk {k = 2 and 3) are constants of integration. Without loss of generality we can set B 3 = 0. At the time wh en we 
match the radiation era solution for 0 to the matter epoch, the argument of the Bessel functions in equation ( 4.36| ) 
is large. For x 1 we have Yi{x) <C Ji{x) so setting i ?3 = 0 is consistent at the time of matching. To see whether 
our assumption is correct we note that the radiation era starts at a time a lot smaller that when it ends. Therefore 
the matching at the other end might be troublesome because for x <C 1 we have Fi(x) :#> Ji(x). To see that we do 
not run into problems we first let x = 4 (a 2 x)^it + l?i)^ so that (redefining Bk , k = 2 and 3, by multiplication by 
an appropriate constant) we get 


0 ( 2 ;) = 


B, 


Y,ix). 


(4.37) 


If we let 0* and 0' be the values of 0(x*) and ^(a^*) at the point of matching (call it x*) to the epoch previous to 
the radiation era, we have 


B3 _ i’*J2{x*) + i’iJijx^) 

B2 0*F'2(x*) + 0(,Tl(x*) 

Therefore even if the matching is done at x* <SC 1 we have 


B 2 


0 {xl) < 1 


at worst. Hence we are justified in taking B 3 = 0 throughout the radiation epoch. Finally we obtain 

B 2 


0(0 = 


(t + Bi )2 


Y Ji[4{a2x)^ {t + Bi 


(4.38) 


(4.39) 


(4.40) 
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E. Epoch of e'^'e Annihilation 


The last era we consider is the one when the electron-positron pairs annhilate. During the early evolution of the 
universe, as we have mentioned, the cosmological fluid can be approximated as a perfect fluid with a radiation equation 
of state. Consider some species i of particles and antiparticles with mass When the temperature of the universe 
passes through the threshold kT ^ rriiC^, the pairs of i particles and antiparticles annihilate, creating important 
though small deviations from the radiation equation of state. These deviations act as a source for our scalar field, 
giving it an extra kick. In particular, we are interested in the thershold for e“'"e“ annihilation, T ~ 5 x 10® K, because 
it is the only annihilation that occurs after the neutron-proton ratio has frozen out of equilibrium, occuring around 
the same time that we have synthesis of primordial elements. 

We will start the e+e" epoch at the time of neutrino decoupling, T ~ 3 x 10^® K. Around such time the universe 
can be safely described as radiation including the behavior of the e“'"e“ pairs because they are still relativistic. We 
will terminate it when the temperature drops to T ~ 3 x 10®. At this point we can once again approximate the 
cosmological fluid as radiation since, even though the leftover electrons are not relativistic their density and pressure 
are very small. Therefore we match with the radiation era solution at this time. 

As derived in Appendix ^ we have the following densities and pressures as functions of temperature: 


1511’ 

(4.41) 

Pi(r) - '’’f ’ 

(4.42) 

Pv{Tv) = -^P-y{Ti^) 

(4.43) 

pAt,) - 

(4.44) 

MT) =1 X - (f 

TT^h Trie kl zrrie kl 

(4.45) 

2mt ,,kT,^^,mes , A:T , , ,2me,, 

PAT) =1 ZPAT) + X )]. 

(4.46) 


Here Ti is Planck’s constant divided by 2tt, k is Boltzmann’s constant, and nie is the electron rest mass. The 
subscripts 7 , v, and e refer to photons, neutrino-antineutrino pairs, and electron-positron pairs respectively. The 
functions Ki{i = \ and 2) are hyperbolic Bessel functions. 

In Appendix ^ we obtain the following differential equations for the scale factor and time as functions of temperature: 


dT\ T 


= 0 


dt 

dT 


dT 


In 


T 

Pi{T)+Pi{T) 


^247T[p,iT) + p,{T)] ’ 


(4.47) 

(4.48) 


where pi(r) = pe{T) +p^ {T ) , pi{T) = pe{T) + p^{T). 

Writing equation (4.16|) for ip using the temperature as the independent variable, we have 




47ra2/2(T)V' = 


(4.49) 


where /i(T) and f 2 {T) are defined by: 


h{T) = 



f2{T) = 


Pi - 3pi 



(4.50) 

(4.51) 


The Ip equation was integrated numerically using a fourth order Runge-Kutta method with adaptive stepsize 
obtained from [^. The initial conditions for the integration were obtained by matching to the radiation epoch 
solution at the lower temperature (T ~ 3 x 10®). There was no need to solve for the scale factor in the e+e" era 
because the speed-up factor depends only on the scalar field and everywhere the scale factor appears in the differential 
equation for ip, it can be replaced by functions of the densities and temperature. 
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F. Nucleosynthesis 


The production of light elements in the early universe depends on the baryon density, the number of light (to ^ 1 
MeV) neutrinos the neutrino lifetime and the speed-up factor in such a way that upper limits on the mass 
fractions X(^He) and -|- ^He) provide an upper limit on To be more specific, the primordial ^He mass 

fraction, X("^He), is given by [ pTpO|] : 

^(^He) = 0.228 -f 0.010 ln(r?io) -t- 0.012(iV,. - 3) + 

0.185 + 0.3271og^„ , (4.52) 


with 7710 = 10^°n(baryon)/n(photon). For the primordial mass fraction of -|- ^He we adopt the conservative upper 
limit X(^H -I- ^He) < 1.7 x 10“"^. This limit implies the inequality 


ln( 7 ;io) > ln(2.60) 0.085(iV,. - 3) + 3.71ogC„ . 

Combining equations ( [4.521) and ( 4.53 ), we find 

X(‘‘He)oi3g^ > ^(^He) > 0.2375 -f 0.013(iVi, - 3) 


0.185 


Ti, — 889.8 
889.8 


0.364 log^n. 


(4.53) 


(4.54) 


We use the values Ni, = 3 and < 885 sec., which are consistent with experimental constraints. We take the 
reasonable limit X(^He) < 0.250, which implies 


InCn < 0.0797. 


(4.55) 


G. Flat and Open Universes 


We now turn our attention to a flat universe (ze ro sp atial curvature, k = 0). For each of these universes we 
calculated the speed-up factor, which from equation (4.17) we know depends only on the value of the scalar held at 
nucleosynthesis. 

In a flat universe, the present epoch can be approximated very well as a matter dominated era. So our solutions 
for the scale factor R{t) and the auxiliary scalar held ip are the ones given in subsection IV C. To evaluate some of 
the constants of integration we impose the appropriate boundary conditions: 


R 

R 


t=to = Ho = 75km/sec/Mpc 
ip{to) = 1 . 


(4.56) 

(4.57) 


These conditions determine two of the constants {k = 2 and 3) in terms of the remaining one. 

Moving back in time we arrive at the radiatio n do minated epoch, which we analyzed in subsection IV D| . We 
can take our radiation era scale factor, equation (4.33), to describe the behavior of the universe up to much higher 
densities, i.e., Bi = 0 , because deviations from radiation dominance (such as annihilation epochs) are short lived and 
only contribute very small corrections. Once again to determine the constants of integration we do the approproate 
matching at the matter epoch radiation era interface, i.e., continuity of i?, ip, and ip. 

Going farther into the past we get to the era of e+e“ annihilation. We have the auxiliary held ip as a function 
of 02 from the solution in the radiation era. Therefore for different values of 02 , we have different v alues of ip at 
nucleosynthesis after the appropriate matching with the radiation era solution. Combining equa tion (4.17), which 
expresses ln^„ in terms of ip{tn), oi, and 02 , with the nucleosynthesis constraint equation ( 4.55| ) for the speed-up 
factor we obtain allowed values of |oi|, for a given 02 which we plot in Figure]^ along with t he so lar system constraint 
(dotted line). The allowed ranges are the regions under the curves. We see from equation ( 4.17|) that for a given 02 , 
|ai| oc 

The hrst thing to note is that throughout most of the range the constraint from nucleosynthesis is more 
restrictive than the solar system constraint. One thing that catches the eye are the peaks at log(a 2 ) — 
—0.55, 1.74, 2.14, and 2.43. Although they appear finite in the graph, they are infinite; meaning that at the 
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peaks ai could have any value (no constraint at all). This is so because, for these values of 02 , the speed-up factor at 
nucleosynthesis is unity (the expansion rate at nucleosynthesis is the same in GR and in ST gravity with these values 
of 02 ). Hence at the peaks the abundances of primordial light elements do not place any restrictions on oi. There 
will be very small ranges around the peaks where the nucleosynthesis constraint is less restrictive the solar system 
one (the first peak is in the region where the solar system constraint is more restrictive anyway). There are values 
of 02 that produce zeros of oi. If there are any deviations (no matter how small) from GR at these values in the 
present, the bounds on primordial light elements will be violated. Most values of 02 are consistent with an upper 
limit of I Oil ~ 10“^ which implies an upper limit on the PPN parameter 7 of |1 — 7 I ~ 10“^°. 

The analysis of the open universe (negative spatial curvature, k = —1) is very similar to that of the flat universe. 
The difference between the two is that the present epoch in the open universe is a curvature dominated one, which 
was analyzed in subsection IV B. The epochs previous to the present one are matter dominated, radiation dominated, 
and a e“'"e“ annihilation era, in that order. We take our open universe to have 0.1 of the critical density, i.e., Go = 0.1. 
For the open universe we plot in Figure ^ the ai vs 02 constraints. The same discussion as for the flat case applies 
here. We have the upper limits |ai| ^ 10“^ and |1 — 7 I ^ 10“^. The main difference with the flat universe results is 
that the limits on oi are 1 to 2 orders of magnitude higher, but still more restrictive than solar system limits over 
most of the range of 02 . 


V. CONCLUSION 

As mentioned in the Introduction, there is motivation to study scalar-tensor theories of gravity from the theoretical 
point of view as well as a means for comparison with experimental gravity. 

After noting the difficulty posed by the arbitrariness in the scalar coupling function a{ip), we noted that its slope 
must be small at present in order to satisfy solar system limits on PPN parameters. In order to remove the freedom in 
the choice of a{(p), we used perturbation theory about the value of the the scalar field today As long as (po is close 
to a value that minimizes a((p), the universe will have general relativistic behavior at present and peturbation theory 
will be equivalent to perturbing about GR, having the slope of a{ip) as the perturbation parameter. In order to find 
nucleosynthesis constraints on oi and 02 we used first order perturbation theory because, if higher order correction 
are impor tan t, they would lead to violation of present bounds on primordial light elements. Moreover, because of 
equation ( [l.3| ) for the speed-up factor ^ ~ 1 indictes that (pn — Po, we can neglect the back reaction on the metric. 

We analyzed the scalar field evolution in flat (Go = 0)and open universes (Go = O.I) to obtain the speed-up factor 
at nucleosynthesis as a function of ai and 02 . Using the constraint that the abundances of light elements place on 
allowed values of ai as a function of 02 were obtained. In both universes the constraint was more restrictive (for 
02 > 1/3) than the one from solar system experiments, the most restrictive one coming from a flat universe. This 
means that if ST gravity is the true theory of gravity, it must be practically indistinguishable from GR through a 
significant period of the universe’s history at least for these cosmologies. 
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APPENDIX A: THERMODYNAMICS OF THE EARLY UNIVERSE 


The purpose of the present appendix is to review and rederive some of the well known results on the 

densities and pressures of different species as functions of temperature. Differential equations for the temperature 
and the expansion factor are also obtained. 

In order to calculate thermodynamic quantities for the different species we need the number density of particles in 
a energy range dE. This quantity is given by 


n{E)dE 


gs PEdE 

exp(^^^) ± 1 ’ 


(AI) 


where gs is the number of spin states, p is the chemical potential, P is the momentum, E is the energy, and the 
-|-(—) sign corresponds to fermions(bosons). For photons and for neutrinos (assuming no neutrino degeneracy) we 
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have fi = 0. For electrons and positrons we assume that the difference between their numbers is small compared to 
their sum, so that we can set fi = 0 for them too (this approximation breaks down when the density and pressure of 
electrons and positrons is negligible). 

The densities and pressures of the different species are given by 


P = 



En{E)dE 


P = 



p 2 

—n{E)dE, 


(A2) 

(A3) 


where m is the mass of the particle species. 

For photons and neutrino-antineutrino pairs (assuming three neutrino species) we obtain 


■K^ikT)* 

(A4) 

P'r 

P, = Y 

(A5) 

Pf ~ '^Plid'v) 

(A 6 ) 

II 

(A7) 


where the subscripts 7 and v refer to photons and neutrino-antineutrino pairs respectively. For the combined densities 
of electrons and positrons we have 


Pe = 



(AS) 

(A9) 


where {i 


1 and 2) are hyperbolic Bessel functions. The following will prove useful in some of our next derivations: 


dpa _ Pa +Pa 

dT ~ T 


(AlO) 


with 0 = 7 , v, or e. 

In the epoch of e“'"e“ annihilation the neutrinos have just deocoupled, hence they do not participate in the thermo¬ 
dynamic equilibrium, but are separetely conserved. The very small baryon contribution is also separetely conserved 
because at these times they have long been decoupled from the equlibriu m. P hotons and e“'"e“ pairs interact and 
are in thermodynamic equilibrium. We have the conservation equation ( [4.13[ )(which is valid to first order in the 
perturbation expansion of the present paper): 


d(pi?3) d{R^) 

dT dT 


(All) 


For neutrinos and for baryons (which we take to be pressureless) the conservation equation leads to 


d{p,R^) 

dT 

djpbR^) 

dT 


(A12) 

(A13) 


with the subscripts v and b refering to neutrinos and baryons respectevely. For the interacting species, e+e pairs 
and photons, we define the following densities and pressures: 


(A14) 
(A15) 

Substituting these last expressions into the conservation equation and using the equality 


Pi — Pe P7 
Pl=Pe=P~f 
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we obtain 


dpi 

d{pi +pi) 

dpi _ 

d{Pi +Pi) 

Pi +Pi 

dT 

dT 

dT 

dT 

T 


3 dR 

^ In/ 

r J, X 



RdT 


,Pl +Pl) ' 



The last equation can be written in the much more compact form 


f R^jPi +Pi] \ 
dT\ T ) ' 


(A16) 


(A17) 


(A18) 


If we combine this result with the field equation (4.12) for the scale factor we end up with 

T 


dt dT 


Pi+Pi 


where 


is the total density. 


dT ^^/24^^p 


p — pi ~\~ Pv Pb — pi Pv 


(A19) 


(A20) 
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FIG. 1. log(|ai|) vs. log(a 2 ) for a fc = 0 flo = 1 universe. Dashed curve is the solar system constraint and the solid curve 
is the nucleosynthesis constraint 

FIG. 2. Iog(|ai|) vs. Iog(a 2 ) for a. k = —1 Do = 0.1 universe. Dashed curve is the solar system constraint and the solid 
curve is the nucleosynthesis constraint 
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